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Abstract 

We study a model of a nonlinear oscillator with a random frequency and derive 
the asymptotic behavior of the probability distribution function when the noise 
is white. In the small damping limit, we show that the physical observables grow 
algebraically with time before the dissipative time scale is reached, and calculate 
the associated anomalous diffusion exponents. In the case of colored noise, with 
a nonzero but arbitrarily small correlation time, the characteristic exponents are 
modified. We determine their values thanks to a self-consistent Ansatz. 
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1 Introduction 

A simple model of a nonlinear stochastic system is obtained by including a 
non-quadratic confining potential in the classical Langevin equation [1,2]. In 
this work, we address the problem of a nonlinear oscillator with random lin- 
ear frequency. Due to the continuous injection of energy by the random force, 
observables such as the oscillator's mechanical energy, position and velocity, 
grow algebraically with time in the absence of dissipation. In Section 2, we 
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calculate the associated growth exponents in the case of Gaussian white noise, 
derive the long time asymptotic behavior of the probability distribution func- 
tion (P.D.F.) in phase space and match it with the (non-Gibbsean) stationary 
P.D.F. in presence of small dissipation. In Section 3, we prove that if the noise 
has a non-vanishing correlation time, the dynamics of the nonlinear oscilla- 
tor is drastically modified: colored noise yields anomalous diffusion exponents 
equal to half the values found for white noise. 



2 White multiplicative noise 



We consider a nonlinear oscillator of amplitude x(t), submitted to parametric 
noise and confined by a potential U(x). For \x\ — > oo, a suitable rescaling 

In . 

permits us to write U ~ t_ with n > 1. Thus, for large amplitudes, the 
stochastic differential equation that governs the system becomes 

|U(f)+ 7 jU (t) + xtf) 2 "- 1 = fflx , with (t{t)£{lf))=V6(t-l!),(i) 



where 7 is the damping rate and £ (t) is a Gaussian white noise of zero mean- 
value and of autocorrelation D. This equation is interpreted according to the 
rules of Stratonovich Calculus [2]. For a quadratic potential (i.e. n = 1) and 
for small 7, the energy of the oscillator grows exponentially with time [3]. The 
effect of nonlinear restoring forces must be taken into account to avoid this 
exponential amplification. 

We shall analyse the motion of the nonlinear stochastic oscillator following 
the method explained in [4,5]. Defining the energy and the angle variables as 

x/E 1 / 2 " 

£ = -i 2 + — x 2n , and = / . , 2 

2 2n (2n) 1 /2« J AfTlp 



we transform the coordinates in phase space from position and velocity to 
energy and angle. Introducing an auxiliary variable Q, with Q = (2n)^~ , 
we rewrite Eq. (1) as 
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where the hyperelliptic function S n satisfies S n ((f) = x/ E l ^ 2n \ These equa- 
tions are rigorous and have been derived without any hypothesis on the noise 
term. In the small damping limit 7 — > 0, we deduce from Eqs. (3,4) that Q ~ t 
and <fi ~ t 2 , i.e. is a fast variable as compared to Q. Averaging the dynamics 
over yields the statistical equipartition identities [4]: 

n 4- 1 

(E) = ^-(x 2 ) and (^) = (i 2 ). (5) 



This averaging procedure [6,7] allows us to derive a closed equation for the 
stochastic evolution of the slow variable Q. We start with the Fokker-Planck 
equation for the P.D.F. Pt(f2, 0) associated with the system (3, 4), and average 
it under the hypothesis that the probability density becomes uniform in <fi when 
t — > 00. Hence, the reduced probability density Pti^l) satisfies 



This averaged Fokker-Planck equation is associated with the following effective 
Langevin dynamics for the variable Q, 

n — 1 \l n + l 



where the effective Gaussian white noise £ has an amplitude V given by 
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i) 2 r £ r 



3 \ 7-1 / 3n+l 



V = V (2n)i ±— )H, 5-^-4. (8) 

" + 1 r 1 r w 



If the dissipation rate 7 is taken to be zero, the solution of Eq. (6) is [4]: 

n + l 

1 n-1 /(2n)^£^\^ f (271)^^1^ 



From this P.D.F. , we derive analytical expressions for statistical averages of 
physical observables [4]. The leading scaling behaviour is given by 

n 1 n 

E ~ (Vt)~ , x ~ (£>i) ^"-d , ^^(Pt) 5 ^ 1 ). (10) 
These formulae are in excellent agreement with numerical simulations [4]. 

n 

We show in Fig. 1 that the averaged energy scales like t™- 1 as predicted. 
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Fig. 1. Scaling behavior of {E(t)) for white, multiplicative noise. Inset: the ratio 
(E) /£ ra /( n_1 ) is plotted vs. time. Eq. (1) is integrated numerically for 7 = 0, T> = 1, 
with a timestep St, and averaged over 10 4 realizations for n = 2, St = 5 10~ 4 ; n = 3, 
St = 51(T 4 ; n = 4,St= 1(T 4 . 



The amplitude x undergoes an anomalous diffusion with an exponent that 
diverges as n — > 1, which is consistent with the exponential growth of the 
linear oscillator. 

In the presence of dissipation, the stationary solution of the Fokker-Planck 
equation (6) is 



n + 1 

1 V 2(n-l)j V 7(n-l)^ / I 7(n-l) ^ J 



This measure is not the canonical Boltzmann-Gibbs distribution. The crossover 
time t c from the asymptotic distribution function for the undamped oscillator 
(9) to the stationary measure in presence of damping (11) is t c — ^^j- • Thus, 
Eq. (9), although derived for a non-dissipative system, is physically relevant: 
it provides an intermediate time asymptotic behavior [8] for the P.D.F. before 
the stationary state is reached; the energy of the system is distributed accord- 
ing to the time-dependent P.D.F. (9) for times t such that 1 <^ t <^ t c and 
then follows the stationary law (11) when t 3> t c . 
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Fig. 2. Scaling behavior of (E(t)) for colored Ornstein-Uhlenbeck multiplicative 
noise. Inset: the ratio (E^/i™/^™ -1 )) is plotted vs. time. Eqs. (12)-(13) are integrated 
numerically for V = 1, r = 1, with a timestep 5t, and averaged over 10 3 realizations 
for n = 2,5t = 1(T 5 ; n = 3, St = 5 1(T 6 ; n = 4, St = 2 1(T 6 . 

3 Colored multiplicative noise 



We now consider a system without damping and submitted to a multiplica- 
tive noise with non-zero correlation time r. In [4], we infered from dimensional 
analysis that the white noise scalings are modified. Here, we shall prove this 
conjecture by calculating the anomalous exponents from a self-consistent ap- 
proach. 

We study the equation: 

^ 2 x(t)+x(tr n -^ V (t)x(t), (12) 

where r)(t) is a colored, Ornstein-Uhlenbeck, Gaussian noise, obeying 

m = -I, (f) + l m . (is) 

at t t 

When t, t' 3> r, we find: 

( V (t)) = and (v(tHt')) = ^e-^. (14) 

As stated earlier, Eqs. (3) and (4) remain valid if £(t) is replaced by rj(t). We 
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assume that, in the long time limit, f2 grows algebraically with time with a 
scaling exponent a, i.e., Q ~ t a . We deduce from Eq. (4) that <fi ~ t u with 
v — a + 1. Substituting this scaling of in Eq. (3), we obtain (leaving aside 
all proportionality constants): 

t 

Sl~J dzSnK)^)^). (15) 


The asymptotic behavior of this expression can be estimated as in [5] by 
discretizing time. This simplification does not alter the critical exponents but 
modifies only the prefactors. We thus replace the colored Gaussian noise rj 
by a discrete random variable e& which takes the values ±^Jv/(2t) randomly 
during the time interval [kr , (k + l)r] . From Eq. (15), we deduce that 

t/r I v 2 t/r 1 

<^ 2 >~E / dz &(z*) £K) ~ £ - 1 3 - 2 ^ , (i6) 

fc=o V £ / k=i (kr) 

where the integral is evaluated by integrating by parts and retaining only the 
leading terms [5]. Assuming that the variable Q is not multifractal, we also 
have fi 2 ~ t 2a . We then deduce that 3 — 2v = 2a. Since v — a + 1, we obtain 
a = 1/4 and ^ = 5/4. Finally, the scaling exponents for colored noise are 

n 1 n 

£^£2(n-l) ; X^t 4 ^" 1 ), X^t 4 ^- 1 ). (17) 



These predictions are confirmed by numerical simulations (see in Fig. 2 the 
scaling behavior of the averaged energy). 

We observe that the exponents are halved when the noise is colored. Further, 
we expect that the colored noise scaling will be observed for large enough time 
as soon as the correlation time r is nonzero. This result may seem paradox- 
ical at first sight because it is often wrongly thought that "in the long time 
limit, colored noise appears white". In fact the period T of the underlying 
deterministic oscillator must be compared with the correlation time r of the 
noise. If r <C T, the noise is uncorrelated over a period of the determinis- 
tic oscillator and acts as if it were white: the scalings found in Eq. (10) are 
satisfied. But, for T < r, the noise remains coherent over a period and its 
effect is perceptible only when it is considered over a large number of periods 
and hence the diffusion slows down: the scalings of Eq. (17) now apply. The 
crossover between the two regimes is observed when T ~ r. The deterministic 

n— 1 

period T decreases with the energy as T ~ E~^~ [4]. From Eq.(10), we have 
_E _2 k- ~ (Dt) -1 / 2 . Therefore, the crossover time scales as t c ~ (Vr 2 )^ 1 . As 
shown in Figure 3, numerical simulations confirm this scenario. Such a crucial 
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Fig. 3. Scaling behavior of (E(t)) for an oscillator in a quartic potential (n = 3) 
submitted to Ornstein-Uhlenbeck multiplicative noise with various values of the 
correlation time. We find a crossover between the white noise behavior observed at 
short times, (E(t)) oc t 2 , and the colored noise behavior recovered at large enough 
times, (E(t)) oc t . 

relevance of the correlation time of the noise was also found in the study of 
the bifurcation threshold of the Duffing oscillator with multiplicative noise [9]. 



4 Conclusion 



A particle trapped in a confining potential and subject to multiplicative noise 
undergoes anomalous diffusion before the dissipative time scale is reached. The 
critical exponents associated with this behaviour are calculated exactly in the 
case of Gaussian white noise and an analytical expression for the intermediate 
time asymptotics of the P.D.F. in phase space is derived. 

Our approach relies on the integrability of the underlying deterministic system: 
we derive exact stochastic equations in energy-angle variables and, after aver- 
aging out the phase variations, a projected dynamics for the energy variable 
is obtained. Our analytical predictions agree perfectly with numerical simula- 
tions. For colored multiplicative noise, the anomalous diffusion exponents have 
been calculated in a self-consistent manner. We find that these exponents are 
halved in presence of time correlations, however small the correlation time r 
may be. 

We observe a crossover from the white noise regime at short times to the col- 
ored noise scalings at long times. The associated crossover time scales like r~ 2 
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and corresponds to the regime when the period of the associated deterministic 
oscillator is of the order of r. In the colored noise case, however, the averag- 
ing technique does not lead to conclusive results because the noise itself is 
averaged out to the leading order. A precise calculation in the case of colored 
Gaussian noise still remains to be done. 

References 

[1] Selected papers on Noise and Stochastic Processes, edited by N. Wax (Dover, 
New York, 1954) 

[2] N.G. van Kampen, Stochastic Processes in Physics and Chemistry (North 
Holland, Amsterdam, 1992) 

[3] D. Hansel and J.F. Luciani, J. Stat. Phys. 54 (1989) 971 

[4] K. Mallick and P. Marcq, Phys. Rev. E 66 (2002) 041113 

[5] K. Mallick and P. Marcq, to appear in Eur. Phys. J. B, preprint 
\cond-mat/0212330\ 

[6] R.L. Stratonovich, Topics on the Theory of Random Noise (Gordon and Breach, 
New York, 1963 (Vol. 1), 1967 (Vol. 2)) 

[7] V. Seshadri, B. J. West and K. Lindenberg, Physica 107A (1981) 219 

[8] G. I. Barenblatt, Similarity, self- similarity and intermediate asymptotics 
(Cambridge University Press, Cambridge, 1996) 

[9] R. Graham and A. Schenzle, Phys. Rev. A 26 (1982) 1676 



8 



